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A  Simulation  of  Weather  in  3D  Space 


I.  INTRODUCTION 

This  paper  describes  a  model  that  stochastically  simulates  a  single-instant 
set  of  values  of  a  weather  element,  both  in  the  horizontal  and  in  the  vertical.  Using 
this  model,  it  is  possible  to  generate  enough  sj-noptic  conditions  to  collectively 
provide  estimates  of  the  probability  distribution  of  weather  events  or  conditions 
in  3U  space. 

In  terms  of  the  element  cloud  cover,  it  is  the  goal  of  this  paper  to  produce  a 
simulated  cloud  cover,  stochastically,  that  will  resemble  natural  clouds,  both  in 
their  horizontal  extent  and  in  their  separation,  to  include  vertical  build-up.  Even¬ 
tually  this  will  allow  us  to  provide  probability  distributions  of  the  amount  of  cloud 
cover  both  in  the  horizontal  and  the  vertical  and  of  the  distances  separating  clouds. 

Stochastic  simulation  is  undertaken  in  the  belief  that  alternative  depictions  are 
not  readilv  available.  Theor'etically  an  abundance  of  observations  in  the  vertical  in 
a  small  anea,  in  rapid  succession,  could  provide  the  kind  of  areal  climatology  that 
■VC  seek.  I’ractiiaillv  speaking,  such  an  abundance  would  confront  us  with  the  pro¬ 
hibitive  task  of  mapping  and  recording  events  as  a  history  of  the  weather,  prior  to 
Its  -.imi  m a  ri/ation  in  statistical  format.  This  is  the  difficulty  that  we  wish  to  over- 
eon.e  bs  modeling.  W C  want  to  he  able  to  simulate  a  sequence  of  changes  in  the 
A.  itln  r  I  iivermg  a  guen  spare  whenever  we  choose,  and  we  want,  ultimately,  to 
pro.ide  11-,  ihh  st  itisties  ineiiiding  correlation  coefficients  (cc)  and  probabilities  of 
Ir,  temal  ,  ,  r. ,  r  1/,  .•■either  conilitions  such  as  cloud  cover. 

1  I I  1  . 1  "  1  I , , '  p I itil  1 1  1 1 1  on  4  •  le  tobe  r  I  '•84  t 
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With  respect  to  time,  most  observations  might  be  considered  as  instantaneous, 
observed  at  regular  intervals.  But  there  are  also  time-integrated  observations, 
like  24-hour  totals  of  rainfall.  There  are  observations  of  the  average  tempera¬ 
ture.  and  maximum  and  minimum  temperatures  in  a  24-hour  period.  All  such 
observations  contribute  to  the  climatological  character  of  a  location  as  it  relates 
the  duration  of  events.  Is  there  a  similar  cMmatic  record  of  areal  or  spatial 
coverage  of  climate? 

In  most  previous  studies,  correlation  coefficients  (cc)  have  been  the  principal 

object  of  investigation,  with  notable  exceptions.  Schreiner  and  Riedel  collected 

2  2 

rainfall  data  integrated  over  areas  ranging  in  size  from  26  km  to  2600  km  ,  and 

found  frequencies  of  extremes  of  rainfall  as  a  function  of  area  of  the  ground  covered. 

2  3  4 

Other  studies,  by  Court,  Briggs,  and  Roberts,  have  linked  the  rain  in  a  circular 
area  to  the  central  single -station  amounts  by  idealized  models.  In  general,  how¬ 
ever.  efforts  to  describe  areal  coverage  have  been  few  and  models  of  areal  cover¬ 
age  have  been  rare. 

Radar  has  been  a  valuable  asset  in  recent  years  because  of  the  panoramic  view 
that  it  provides  of  precipitation  within  a  circular  area.  Visual  satellite  imagery, 
recording  as  it  does  the  areal  cloud  cover,  eventually  should  help  to  provide  a 
21)  clcjud  id iniatology.  Satellite  IK  gives  the  tops  but  not  the  internal  structure. 
Satellites  and  radar,  however,  .are  recent  toots.  Kxisting  climatic  records  gener¬ 
ally  contain  single -instant  obsei'vations  at  single-point  weather  stations.  A 
re.'oried  rainfall  amount  is  that  at  a  r.ain  gauge.  Temperature  is  that  in  a  ther- 
nama'ter  shelter,  and  so  on, 

A  word  shout  1  be  sal  1  about  numerical  weather  prediction  (NWP).  It  would 
-e.'m  that  Ntt  1’  offers  .in  appi-oach  to  our  problem,  to  generate  a  sequence  of 
instantaneous  hon/ontal  pictures  at  rt.-gular  time  intervals.  This  may  be  done  in 
die  future.  At  this  point  in  time,  however,  stationarity  of  NWP  models  in  a 
-tatistii  al  sense,  is  considered  unproven  and  doubtful.  The  use  of  NWP  to  esti¬ 
mate  .  limatic  nrobalu  1  Hies  is  expensive  and  time  consuming,  yet  deficient  in 
supiilving  climatologies  of  the  items  of  interest,  including  cloud  cover,  ceiling, 
an  I  \  isihility,  IVy  and  large,  we  must  depend  on  statistical  models. 

1 . 


Sctireiner,  1..  f.  .  and  Riedel.  J.l'.  (I't78)  Probable  Maximum  Precipitation 
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2.  EQUIVALENT  NORMAL  DEVIATE  (END) 

Before  introducing  the  model,  we  need  to  examine  the  concept  of  the  equivalent 
normal  deviate  (END).  Single -point  frequency  distributions,  for  a  variety  of 
weather  elements,  have  been  studied  in  a  variety  of  ways.  Generally  they  do  not 
have  a  normal  distribution.  But  there  always  is  a  value  (END)  of  the  normal  dis¬ 
tribution  that  has  the  same  cumulative  probability  as  any  element  value.  For 
example.  Figure  1  shows  the  cumulative  relative  frequency  of  visibility,  specific¬ 
ally  for  the  January  noontime  visibility  at  the  Bedford,  Mass,  airport.  Somerville 

5 

et  al  have  fitted  the  visibility  data  to  the  Weibull  distribution,  as  follows: 

F(V)  =  1  -  exp(-aV^)  (1) 


where,  for  Bedford,  January  noontime. 


a  =  0.  06906 

,i  =  0.  8186  for  V  in  miles  . 


Figure  1.  Cumulative  Relative 
Frequencv  of  the  January  .Noontime 
X’isibility  at  the  Bedford,  Mass. 
■Airport.  The  x's  mark  the 
frequencies  as  summarized  from 
20  vrs  of  data.  Thev  show,  for 
example,  a  2  percent  frequency 
for  visibilitv  less  than  1/4  mile, 

7  percent  for  visibilitv  less  than 
1  mile,  and  38  percent  for  visibility 
less  than  10  miles,  that  is, 

62  percent  frequencv  of  visibility 
equal  to,  or  greater  than,  10  miles 


» 


i 


» 
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5.  Somerville,  F.  N.  ,  Bean,  S.  J.,  and  Falls,  S.  (1970)  Some  Models  for  \'isibilitv. 
Scientific  Report  No.  3.  Contract  No.  F 19628 -7  7  -C -0080 ,  Universitv  of 
Central  Florida,  Orlando,  Florida,  .A  FGL -TR -7t< -0  144,  AD  A073490. 
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The  vuniulativi'  relative  ('I'equeney,  i-'(V).  ranges  from  0  to  1.  The  END-scale  of 
i'(\,  I  is  [ilntted  as  v  on  the  graph  alongside  the  probability  scale,  and  ranges  from 
-  ^  to  ■  .  Thus  we  obtain  a  transformation  of  visibility  (V)  into  its  END  (y)  via 

l''(\  1  or  vice  versa. 

Stochastic  models  derived  using  the  END'S  instead  of  the  original  meteorological 
\  ariables  tend  to  be  more  robust  and  applicable  to  a  greater  range  of  situations. 

Tile  results  arc  transformed  to  the  appropriate  weather  variables,  when  so  required. 

If  we  are  given  the  cumulative  relative  frequency,  F(x)  of  any  variable  (x),  it 
is  defined  also  as  the  cumulative  relative  frequency  of  the  END  (y).  That  is, 

K(v)  “  K(x)  . 

g 

To  calculate  y  given  F(v),  the  National  Bureau  of  Standards  provides  the 
afiproximate  ffirmula 


w  he  re 

a  =  :i.:i07  5:i 
o 

aj  -  0.:a7061 
bj  --  0.99229 
b^  ^  0.04481 


■where 

p  ■  F(v)  . 


0.  N  HS  (1984)  Handbook  of  Mathematical  Functions,  .■\pplieri  Matliematical  Senes, 
No.  n.9.  (lovernment  Printing  Offire,  Washington,  D.  C. 


U  hen  the  I’.M)  (v)  is  provided,  the  corresponding  cumulative  frequency,  F(y) 
is  afjproxiri'.ated  bv  an  NBS  formula,  as  follows: 

l'(v')  =  <  -  m  j  2(1  +  Cj  'y'  +  c^y^  +  c^  'yi"^  +  |  ^ 

whe  re 

Cj  -  0.  19  6854 

c,,  -  0.  115194 

c.^  =  0.000344 

c,  ^  0.019527 

4 

/  -  0,  m  =  1  for  y  <  0 

'■  ■  1.  m  -  - 1  for  y  ?  0  . 

■file  cumulative  relative  frequencies  from  the  observed  sample  are  considered 
estimates  of  the  l  umulative  probabilities  of  the  true  distribution.  If  the  operative 
svnibol  for  probability  is  P,  we  assume: 

P(y)  -  Plx)  =  Fix)  =  P(y)  .  (4) 

Another  exaniftle  of  the  transformation  of  a  meteorological  el  n  •'nt  into  its 
PM)  is  shown  in  I'igure  2.  The  MND-scale  corresponding  to  the  Cu...ulative  rela- 
'ive  ir'equencv  (hence  the  probabilitv  estimate)  of  24-hour  single -station  precipita- 
'lon  IS  shown  f(jr  New  Mngland.  during  the  month  of  January. 

J.  Of  SIRED  MODEL  ATTRIBUTES 

1  o  hecin  with,  the  process  of  generating  a  model  of  the  areal  extent  of  the 
•w  iflier-  can  resemble  a  Markov  process,  but  it  cannot  be  one.  In  a  Markov  model 
I  i  n'nre  I'vent  is  defiendent  upon  the  present  state  of  the  weather  and  is  not  addi - 
lonallv  governed  by  events  prior  to  the  present.  In  a  horizontal  picture,  any  one 
ooint  '  At  is  surrounded  by  a  double  infinity  of  points.  There  must  be,  consequently, 
^ome  iinderlyiriL!  restraint  on  the  conditions  surrounding  A,  similar  to.  but  still 
lifferent  from,  a  Markov  process. 


24- HOUR  RAINFALL  (MM) 

FiL’ure  2.  Twenty-four  Hour  Single -Station  Precipitation  for 
New  Fngland,  Based  on  27  Januaries  (1952-1978).  The 
ordinate  is  the  uniform  scale  of  y,  the  END.  The  value  for 
zero  preci pitation  is  y?0.4,  for  10  mm  or  less  it  is 
V  1.4,  and  so  on 

An  effective  model  of  areal  coverage  will  produce  a  simulated  horizontal 
field,  slochasticall V,  that  im.itates  nature  but  without  direct  recourse  to  physical 
laws  or  dvnamics.  .An  example  of  the  kind  of  field  that  the  model  should  simulate 
IS  shown  in  {-'igure  3.  Once  a  model  is  selected,  the  adopted  procedure  can  be 
r  eneated  m.any  times  to  produce  manv  simulated  synoptic  situations,  thus  com - 
iirising  a  large  sample  that  can  be  surveyed  for  spatial  correlations  or  associations 
)f  c')n  iitions,  Therefoi  the  miodel  must  be  fast. 

Idle  coeffii dents  of  correlation  of  values  at  two  points  in  the  space  should 
generall  v  iet  rease  with  increasing  distance  between  the  two  points  in  the  same 
manner  as  nature.  Hertoni  and  Lund*  have  presented  results  for  cc's  of  pressure, 

7.  iiertoni.  If.  .A.,  and  Lund,  I.  A.  (  1964)  U  inter  .Space  Correlations  of 
I’cessure,  l  emperature  and  Density  to  16  km;  .Al-'C  1<L -64  -  1020, 

AD  6  1  1002. 


(  )bric  i'\  <ti<>na  1  info  filiation  on  cloudiness  has  a  reasonably  high  degree  of 
accurac'.  ,  as  long  as  it  is  limited  to  the  determination  of  ceiling,  or  the  U’vel  at 
which  the  skv  is  cuimilatively  6/10  covered.  There  is  more  information,  however, 
in  the  daily  and  hourlv  observations.  In  the  TD1''-14  Airways  Observations,  the 
amounts  of  clouds,  in  lavers,  up  to  four  layers  aloft,  have  been  recorded.  Ob¬ 
servers  must  face  the  dilemma  that  lower  clouds  block  the  view  of  higher  clouds. 

It  is  possible  to  confidently  record  the  amount  of  a  low  cloud.  Hut  the  amount  of 
cloud  in  a  higher  layer  could  be  as  low  as  the  extra  amount  of  sky  covered  by  the 
higher  cloud,  or  could  be  as  high  as  the  total  sky  covered  up  to  that  level. 


Table  2. 
a  Ceiling 

The  A  re 
of  Clouds 

il  Coverage  A  Corresponding  to 

Whose  Base  is  at  11  (Eq.  (25a)] 

Cloud  Height  (11) 

Areal  Coverage  (A) 

(ft) 

(km) 

km^ 

0.  1 

1.  0  (arbitrary) 

500 

0.  15 

2.  3 

1, 000 

0,  30 

9.  1 

o 

o 

o 

0.  91 

83.  3 

5,  000 

1.  52 

232 

10, 000 

3.  05 

925 

20, 000 

6.  10 

3,641 

o 

o 

o 

o 

9.  14 

8,  050  km^ 

The  Hnvironmental  and  Technical  Applications  Center  (ETAC)  of  USAF  made 
the  latter  assumption  to  provide  a  distribution  of  the  subtotal  cloud  cover  for 
12  levels  from  the  surface  to  30,  000  ft  (Table  3).  As  long  as  the  item  of  interest 
is  the  sub-total  cloud  cover,  from  the  ground  up  to  a  given  level,  the  assumptions 
about  the  hidden  clouds  are  not  relevant.  As  the  height  is  increased  the  probability 
of  a  clear  sky  (0/H)  decreases  monotonically,  while  the  probability  of  (8/8)  increases 
If  f  IS  the  symbol  for  fractional  cover  (f  =  0/8  (1/8)  8/8],  the  cumulative  probability 
of  the  fraction,  l’(TT),  that  is,  the  amount  (f),  will  decrease  monotonically  with 
height.  (This  has  been  verified  in  Table  3.) 


VERTICAL  CLOUD  SIMULATION 
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Figure  9.  Another  Stochastically  Generated 
Picture  of  a  V'ertical  Cross-Section,  With  the 
Same  Conditions  as  in  Figure  8,  but  Producing 
Little  Cloudiness 
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Figure  10,  Another  Stochastically  Generated 
V'ertical  Cross-Section.  Which  Depicts  Two 
Cloud  Layers 
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Figure  1 1.  An  Additional  Stochastically  Generated 
V'crtical  Cross-Section 


5.3  Portrayal  of  Mean  Cloudiness— in  Layers  and  at  Levels  Aloft 

As  seen  above  it  is  necessary  that  we  specify  the  climatology  of  cloudiness  in 
selected  layers  or  levels  aloft.  This  section  describes  a  method  to  estimate  the 
vertical  climatology  of  clouds  given  standard  weather  observations. 


VERTICAL  CLOUD  SIMULATION 
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Figure  8.  Sample  of  a  Stochastically  Generated 
Picture  of  a  Vertical  Cross-Section  of  Clouds. 
Chosen:  A  =  256  km  in  the  Horizontal. 
Shrinkage  in  the  vertical:  1/50.  Mean  cloud 
cover  (Pq)  varies  from  0.58  (Yq  =  -0.20)  in  the 
lowest  1.  5-km  layer  to  0.  06  (y^  =  1.  56)  in  the 
highest  1.  5-km  layer.  Cloud  presence  is 
indicated  whenever  the  simulated  END(y)  is 
greater  than 


The  mean  cloud  cover  in  the  first  layer  at  the  location  (Central  Europe  in 
.Januat'v)  is  0.58,  to  give  the  Iv-N'D.  v  -  -0.20.  The  x's.  therefore,  in  Figure  8 
were  located  wherever  v  5  -0.20.  In  the  second  layer,  from  1.5  to  3  km,  the  mean 
cloud  cover  is  0.  55,  to  give  y  =  -0.  13,  and  so  on.  For  the  top  layer,  from  7,  5  km 
to  9.  0  km,  the  mean  cloud  cover  is  only  0.06  corresponding  to  y  S  1.56,  hence 
there  is  infrequent  cloudiness  in  this  layer. 

.-\s  opposed  to  the  overcast  of  Figure  8,  another  picture,  stochastically  pro¬ 
duced  with  the  same  wavelength  and  mean  cloud  frequency,  resulted  in  the  relatively 
clear  sky  of  Figure  9.  Still  another  example  (Figure  10)  produced  two  layers  of 
clouds,  whose  total  cloud  cover  was  0.73,  compared  with  the  mean  total  cloud 
cover  for  Central  Europe  of  0.  7  8.  .Another  example  (Figure  11)  was  produced, 
wherein  the  random  numbers  that  were  used  for  the  horizontal  picture  of  Figure  7 
Were  deliberately  used  for  the  vertical  cross-section  in  Figure  11. 

The  :tbove  examples  (1- igures  7-11)  are  aimed  at  producing  a  picture  of  a 
::loud\  field,  t.'ither  in  the  horizontal  or  in  a  vertical  cross-section.  Alternatively, 

It  i.s  possible  to  imagine  th.at  :i  vehicle  s  line  of  travel  is  plotted,  and  a  determina¬ 
tion  HKide  as  to  where,  ;ilong  it.s  path,  it  penetrates  a  cloud.  So  far,  it  would  have 
to  be  assumed  that  the  vehicle  completes  its  travel  instantaneously,  or  within  a 
reasonablv  short  intei'val  of  time,  to  avoid  considering  time  changes. 


then  the  cloud  boundaries  should  be  found  wherever  the  KNIJ  y  -  -0.26.  Hence,  the 

shading,  which  represents  the  cloud  presence,  has  been  inserted  on  the  figure 

wherever  y  =  -0.  26.  Three  circles  of  radius  27.  8  km  are  centered  at  three 

supposed  observation  points,  where  the  observer's  sky  dome  has  a  radius  of 

2 

27.  8  km  (15  nm)  or  a  floor  space  of  2424  km  .  At  one  of  the  three  areas,  the 
observer  would  see  a  clear  sky,  at  another  a  broken  sky,  and  at  the  third  he  would 
record  an  overcast. 

Table  1.  A  Typical  Set  of  Values  for  hj,  cos  Oj,  cos 
cos  yj.  i  =  1,  12  That  Were  Obtained  by  Random  Process 
and  Used  to  Obtain  Values  of  y(u.  v,  w)  by  Eqs.  (17), 

(16),  (19) 


i 

h. 

1 

COS  a. 

1 

COS 

1 

COS  V  . 

'  1 

1 

0.  3036 

0.  5159 

-0.  0994 

-0.8509 

2 

0.4705 

-0.  2436 

0.  9699 

-0.  0050 

3 

0.  1 159 

-0. 1818 

0.  0045 

0.  9833 

4 

0.  0944 

-0. 1183 

-0.  0069 

-0.  4044 

5 

0.  9345 

0.  2931 

0.  8873 

-0. 3660 

6 

0. 6600 

-0. 1013 

-0.  7891 

0.  6059 

7 

0.  5655 

0,  1346 

-0, 3544 

-0. 9253 

H 

0. 3212 

0.  9982 

-0. 0366 

-0. 0479 

fl 

0. 6899 

-0.  2305 

-0. 6861 

0, 6900 

10 

0.  3618 

0.  1872 

-0. 9390 

-0. 2885 

1  1 

0.  2223 

-0. 8925 

-0. 2022 

0.4031 

12 

0. 6509 

0.  3380 

0.  4760 

-0. 8119 

5.2  Sample  Vertical  Cross-section 

T'l  L’cnc  iMtc  th('  vt  i'tK;:il  picture  of  Figure  8  the  wavelength  was  assumed  to  be 
\  -  256  ki:.,  Hi  101250  kni,  v  ;  0  (for  a  single  cross-section).  In  the  vertical 

the  elo'irl.^  ueni- r  . ted  for  .six  l:iyer.s  of  1.5  km  thickness,  to  make 

=  1 .  5 1  1 .  5  )  '  ■ .  0  r.  I ! .  . 

Ttii-  1..\IJ  -T  f'lund  ..t  26  succes.sivt-  points  horizontally  for  each  of  the 

fj  1  ivtM'.s  iloft.  li  '.  i.s  i.-,,suti,cd  thiit  1.5  km  in  tlie  vertical  was  equivalent  to  75  km 
in  tlie  noi  i/cnt  !.  l';im  v, ,  for  our  purpose.s  in  model  application;  w  =  q.  w'. 

In  tlii.T  '  Xiirplr,  u  511,  t.o  L'l'vc  ',,5'T5)450  km. 


5.1  Sample  Hori^unlal  Simulation 


111  tills  (1  ^  iOO  kni,  I  -  10,  w  -  0,  and  u,  v  were  each  given 

1  i  .s Lie ct' s s i\'t‘  \'alLU‘S; 

u  0l2."'i)  .100  kni 

V  -  0(2,'))  200  km  . 


Figure  7.  Sample  of  a  Stochastically  Generated  Field  of 
Clouds  in  a  Horizontal  Plane.  Chosen:  A  =  300  km, 

Py  =  0.6,  so  that  clouds  are  depicted  wherever 
V  =  -0.  26.  Solid  lines  indicate  constant  ENDs;  hatching 
indicates  area  covered  by  clouds.  Circles  indicate 
three  possible  areas  of  radius  27.8  km  an  observer  might 
see 


Using  a  random  number  generator,  a  total  of  12  X  3  =  36  random  numbers  were 
generated,  and  then  transformed  into  the  numbers  (Table  1)  that  must  characterize 
the  resulting  field  of  END-  or  y-values.  Figure  7  shows  the  results  of  drawing 
isopleths  of  the  END'S  (y)  between  the  (13  X  13  =  169)  values  generated  at  the  169 
grid  points.  If,  as  at  Bedford,  Mass.  ,  in  January,  the  mean  cloud  cover  is  6/10, 


T" 


T 


r','  V' 


11'  this  fii-lti  is  modified  by  addinfi  a  second  wave  formation  to  the  first,  then  a 
pattern  Aill  emerge  that  will  be  repeated  approximately  at  1  or  2  wavelengths.  If, 
!in.\e\er,  a  tliii'd  wave  is  added  its  period  will  tilmost  alw'ays  (depending  on  the 
random  number)  be  an  irrational  fraction  of  the  period  of  the  first  two  waves.  As 
more  waves  are  added  the  prijbability  of  all  the  waves  having  rational  periods  in 
ari.-  ri\eii  lireeij  iti  is  infinitesimal.  Ibis  intplies  tliat.  with  many  wave  formations 
(  ]2  or  ’!!  tre)  the  (itiitern  iti  anv  otie  small  area  is  not  likely  to  be  repeated  or  relateil 
a  anv  w  i  /  to  inotlier  lisltitii  ptitterti. 

5.  APPLICATION  TO  SIMULATION 


■('he  field  elut racte ristic  (A)  is  presumed  known,  and  a  simple  value  (like  10) 
is  assigned  to  the  integer  (1).  The  stochastic  procedure  requires  that  N  (=  12,  say) 
wave  formations  be  generated  stochastically,  which  means  finding  N  sets,  each  of 
J  random  numbtjrs,  to  obtain  values  |Eq.  (15))  for 

h, .  cos  ■}  . ,  A , .  i  =  1 ,  N  . 
i  11 

which,  in  turn,  yield  the  v.alues  |  l.qs.  (H)  and  (!))]  for 


h.,  cf)s  (».,  ttos  cos  i  =  1,  .N  , 

I'ttr  the  point  (It,  V,  w  )  in  spticc  lUq.  (17)] 

I)  ;  u.  cos  <y .  ^  V.  cos  ,i.  +  w.  cos  V  . 

1  i  '  i 

to  give 

X  -  1  HA  (h.  r  U  'A  +1)  i  =  1,  \ 

1  1 

to  \ield  the  M.\U: 


,\ 


(21) 


(22) 


(23) 


I'o  obtain  a  horizontal  |iicturc,  stochastically,  the  distances  (u  and  v)  should 
hi’  t  hosen  at  regular  gihd  points,  while  the  vertical  height  (w)  should  be  kept  fixed. 
Say  at  w  :  0.  To  obtain  a  viTtical  picture,  one  horizontal  dimension  (v,  say)  should 
be  ke[)t  fixed,  while  (u,  w)  are  assigned  regular  gridpoint  distances.  The  result 
will  he  a  ■  foss-sectional  view  in  the  U-\V  plane. 


14 


Suppose,  now,  that  there  are  N  formations.  Then,  at  each  point  (u,  v,  w)  there 
will  be  a  set  of  N  wave  phases,  whose  average 
N 

X  =  1/N.  X]  x.(u,  V,  w)  (18) 

i=  1  ^ 


will  have  a  mean  (1/2)  and  variance  (1/12N).  The  standardized  value  at  any  point, 
y(u,  V,  w),  is  given  by 


y(u,  V.  w) 


(19) 


The  sum  of  N  uniform  variables  is  used  to  approximate  a  normal  distribution;  ^ 
therefore  y  closely  approximates  an  equivalent  normal  deviate  (END). 

There  will  be  a  value  for  y  at  each  and  every  point  in  the  three-dimensional 
space.  Thus,  beginning  with  3N  random  numbers,  it  is  possible  to  generate  a  map, 
or  condition  or  shapshot  picture,  of  END-values  in  the  space  to  resemble  isopleths 
of  temperature,  or  rainfall,  or  cloud  cover,  and  so  on,  at  various  levels  in  the 
atmosphere. 

When,  for  convenience,  N  is  made  equal  to  12,  the  equation  becomes 
12 

v(u,  V,  w)  =  x.(u,  V,  w)  -  6  .  (20) 

1  =  1  '■ 


For  clouc  cover,  the  wavelength  (A)  will  be  several  hundred  kilometers  in 
length,  while  for  pressure  or  tempoi-atuie  !iel  ls  n  will  be  as  much  as  one  thousand 
kilometers  or  more  m  length. 

A  vertical  dimension  of  20  km  should  correspond  to  a  horizontal  dimension  of, 
say,  1000  km.  Th;it  is,  if  we  trv,  for  example,  to  simulate  a  horizontal  field  at 
sea  level,  and  also  the  field  1  km  above  sea  level,  we  might  do  so  by  scaling  to  an 
effective  separation  of  aO  km.  The  actual  equivalence  is  yet  to  bo  determined. 

A  question  arises  about  the  recurrence  of  patterns  in  a  given  map.  In  the 
horizontal,  will  the  picture  produced  in  one  small  area  of  the  field  bo  repeated  in 
another  area  ’  If  the  procedure  stopped  with  one  wave  formation,  then,  as  in 
Figure  5,  there  would  be  a  regular  pattern,  with  wave  phases  repeating  at  distances 
of  1.  2,  - wavelengths. 


9.  Enslein,  K.  ,  Ralston,  A.,  and  Wilf,  H.  (1977'  Statistical  Methods  for  Digital 
Computers,  John  W  ilev  and  Sons,  Inc.,  Ne  .  York. 


Also  since  each  element  is  directly  proportional  to  cosy  ,  we  select  cosy  uniformly; 


-  1  2  cos  >21.  (13) 

A  wave  formation  that  fills  the  3-dimensional  space  will  have  three  degrees 
of  freedom  (two  for  the  angles  and  one  for  h,  the  phase  at  the  origin).  The  phase 
will  be  uniform  between  0  and  1, 

0  2  h  2  1  .  (14) 

Hence,  to  utilize  the  three  degrees  of  freedom  for  the  selection  of  a  wave 
formation,  three  random  numbers  (H^^.  ai'e  selected  to  give  h,  cosy.  A), 

each  with  a  rectangular  distribution,  as  follows: 

h  = 

cos  >  =  2  R  -  1 

> 

A  =  2,;  R^  .  (15) 

U  ith  a  single  wave  formation  there  is  a  determinable  value  of  the  wave  phase  (x) 
at  each  point  (u,  v,  w),  as  follows: 

X  =  I'ltA  (h  +  U/A  +  I)  (16) 

where 

U  =  u.  cos  a  +  v.  cos  ,i  +  w.  cosy  ( 17 ) 

and  where 

cos  (Y ,  cos  ,(  are  found  by  Mqs.  (8)  and  (9)  . 

Now  suppose  that  another  wave  formation  (with  the  same  wavelength  A)  is  pro¬ 
duced  through  the  choice  of  three  other  random  numbers.  Again  there  will  be  a 
determinable  value  of  the  second  wave  phase  at  each  point  (u,  v,  w)  but  the  sum,  or 
!',cr age,  of  the  tw(j  values  will  show  a  more  intricate  structured  pattern  than  those 
of  one  w  ive  formation. 


w 


r 


i'igure  6.  A  3-Diinensional 
View  of  a  Point  (P)  With 
Cartesian  Coordinates  (u,  v.  w), 
at  Distance  D  Prom  the  Origin 
(0,  0,  0),  With  Spherical  Sides 
(a,  >  )  and  Angular  Measure 

A  Prom  the  U-W  Plane 


I’roni  the  spherical  triangles. 


cos  (>  -  sin  7  .  cos  A 


cos  ,i  =  sin  >  .  sin  A  . 


2  ■;  ■> 
cos  a  -  cos“,<  •  (:os“>  =  1 


w‘iii  1;  1!:  rjlics  tha'  a r' ill t ra r '/  ■  iioii’cs  for  sides  ((V,  /3.  >)  arc  effectivel'/  reduced  to 
•  A  O.  I:'  '  an  i  )  are  sold  te  i,  (V  and  (i  can  he  calculated  bv  Pqs.  (8)  and  (0). 

\  sphere  ol  radius  1)  has  .i  surface  area  of  4  .t  U  .  An  clement  of  the  area  (6A) 
of  tins  S[)!iere  is 

'  A  ;  (  D.  1  7  )  ( U.  sin  7  .  "lA  ) 

2 

D“  .  5  ( -  cos  7  )  .  ^)  A  .  (11) 

Ac  ■.■..,nt  each  element  to  be  equally  probable.  Since  each  element  is  directly 
pr  oportional  to  fiA  ,  sede.  r  unil'ornilv; 

'J  ■  »  •  3  . 


(1  =  u.  108  0  ■  V.  sin  0  , 
which  can  h>  negative  on  positive. 

An  integei-  (I  '  10,  sav)  is  added  to  the  sum  (h  •  d'A)  to  assure  a  positive 
\alue  for  x. 


(6) 


1 


l-'igure  ■).  A  I’lan  \  iew  and  (.' ross -Section  of 
a  S.ivvlootli  Wave  I'ormation 


lor  the  :l|)  HS\\  no  del.  coitsider  ativ  point  (I’l  with  the  three  coordinates 
(u’,  v'.  rt').  A  lineCrom  theorigin  toP(ii',  v'.  w' )  is  oriented  wi  th  respetd  to  the  origin 
(0,  1,  0i  bv  angles  ui ,  •  between  the  line  tmd  the  u,  v,  w  axes  respectively  as 

shown  in  1  igiiie  fi.  consider  the  plane  perpendicular  to  the  line  at  1’.  The  per¬ 
pendicular  distance  of  the  plane  fiom  the  origin  is  L).  I'rom  solid  geometry, 
the  equation  of  the  (ilane  is 

Li.costy  +  v.cos.C  ■>  w.eosy  IJ  .  (7) 


Let  A  be  the  angular  measure  from  the  I  -VV  plane  to  the  point  P.  (If  the  W-axis 
were  the  earth's  axis,  A  would  be  the  longitude.) 


Hertoni  anci  Lund  found  the  cc's  for  the  readily  available  elements  of  pressure 
and  temperature,  which  displayed  much  persistence.  The  spatial  correlations  for 
(Phenomena  sui'h  as  precipitation  and  clouds  are  expectedly  less.  The  zero  correla¬ 
tion  distances,  compared  to  those  in  Figure  4,  should  be  much  shorter.  An 
effci  tive  model  will  produce  a  relation  of  correlation  with  distance  that  will 
resemble  those  found  from  the  raw  data. 

Tiic  model  needs  to  be  able  to  give  the  probability  of  joint  events  or  events 
occurring  in  given  fractions  of  an  area  or  space.  For  example,  cloud  cover  ranges 
from  I  dear,  through  partial  cover,  to  overcast. 


4.  THE  THREE-DIMENSIONAL  BOEHM  SAWTOOTH  WAVE  (3D-BSW)  MODEL 

d'he  objei dive  of  the  model  is  to  produce,  stochastically,  on  each  map  or 
occasion,  a  correlated  field  of  values  of  the  RND  similar  to  the  field  in  Figure  3, 
and  to  produce  such  a  depicdion  at  any  level,  or  to  produce  cross-sectional  pictures 
in  the  vertical. 

liegin  with  a  stationary  wave  formation  that  fills  the  space.  Figure  5  shows  a 
[dan  view.  de()i<ding  the  waves  as  having  a  uniform  wavelength  (A).  (The  formation, 
at  this  stag(>  of  the  explanation,  is  limited  to  the  two  dimensions  of  the  figure.  ) 
d'he  leariing  <‘dge  of  the  middle  wave  misses  the  origin  of  the  U,  V -axes  by  the 
rlistam.'o  (ID,  that  is,  the  phase  (h)  at  the  origin  is  h=  H/a  and  0  •?  h  1  1.  The 
orientation  of  the  waves  is  given  by  the  angle  (0).  The  values  (h,  ff)  are  randomly 
chosen,  but  the  wavelength  (A  km)  is  a  parameter  of  the  model  whose  length  (km  or 
nm)  depends  on  the  fluctuation  of  events  in  space. 

The  sawtooth  nature  of  the  waves  is  depicted  by  their  cross-section  (bottom  of 
I'  igure  5).  The  front  edge  of  each  wave  has  a  height  of  0,  the  back  edge  has  a 
height  of  1  and  in-between  points  have  a  uniform  distribution  of  heights  (x)  where 
X  is  a  number  from  0  to  1.  That  is,  the  height  of  the  wave  equals  the  phase.  Thus, 
at  eac  h  point  (u.v)  of  the  space,  there  is  a  wave  height  (or  phase)  (x): 

X  =  IdtA  (h  •  d/A  •  1)  (5) 

where  I  MA  stands  for  the  fractional  part  of  the  bracketed  sum,  and  d  is  the  dis- 

g 

ance  (kill)  of  the  wave -phase  line  from  the  origin,  in  the  two-dimensional  model 
■he  distance  (i|)  is  found  to  be 


'1.  Miircei'.  ('.1-.  .  and  (i ri ngorten,  I.  I.  ( 11)84)  Two-Dimensional  Modeling 

lor  Lineal  and  Areal  Probabilities  of  Weather  Conditions.  A  FCi  L-TK -84 -0  126, 
AT)"YTT7T77T - - - 
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Table  3  has  the  necessary  data  to  find  the  mean  cloud  cover  from  surface  to  a 
given  height  (last  column),  which  increases  monotonically  with  height.  This  kind 
of  information  permits  the  drawing  of  one  curve  in  a  diagram  of  mean  cloud  cover 
versus  height  (Figure  Ik).  But  there  is  still  a  need  for  more  information,  for  the 
cloud  cover  in  a  layer. 


30,000 
Ht 

(tf)  20,000 
10,000 

0  I  2  3  4  5  6  7  8  ,9i0 

mean  cloud  cover 

Figure  12.  The  Wean  Cloud  Cover,  From 
the  Ground  Up  to  Given  Heights,  Fort  Worth, 

Texas,  January,  0300L 

DeBary  and  Woller^*^  obtained  valuable  information  on  the  cloudiness,  in  layers 

aloft,  through  observations  made  from  aircraft  flying  over  Germany.  Their  study 

remains  unique  with  respect  to  the  source  of  data. 

For  a  more  universal  approach  ETAC  has  used  the  data  of  3D-Nephanalysis  as 

gathered  and  processed  by  Global  Weather  Central.  The  latter's  sources  include 

worldwide  surface  observations  and  satellite  data.  ETAC's  data  on  cloud  cover 

■  listributinns  were  submitted  to  us  for  4.S  layers  between  10  levels,  taken  in  pairs 

the  .surface  to  55,000  ft.  Table  4  gives  a  sample  frequency  distribution  of 

'  lou  i  aniount,  p(fl,  f  -  0(0.  1)  1.0  as  estimated  within  the  specified  layer. 

'I'abic-  4  IS  one  of  (14  /  45)  tables,  one  for  each  of  7  locations  (Korea,  north 

mil  a,t)uth;  Central  Europe,  .Northern  Plains  and  Southern  Mountains;  Mid-East, 

Iji  -,1  I  t,  Eastern  Mediterranean  and  Mountain  climates),  for  each  of  two  mid- 

i-ion  nionths  (January  and  July)  and  for  each  of  45  layers  in  the  atmosphere. 

■-mi  ii  frequencv  distribution  of  cloud  amounts  vield  the  mean  cloud  cover  (P  )  bv  the 

o 

estima  te: 

1.0 

P  P'f  ■  f  f  =  0(0.  1)1.0.  (24) 

r  =  0 

10.  IJeHarv  1C.,  and  Moller,  E.  (1963)  The  vertical  distribution  of  clouds, 

J.  .Appl.  Meteorol.  ,  2:806-808. 


Table  4.  Example  of  the  Frequency  Distribution  of 
Cloud  Amount,  p(f).  Within  a  Layer  of  the  Atmos¬ 
phere.  The  example  is  for  Central  Europe,  Northern 
Plains,  January,  layer  from  1000  to  to  3500  ft 
(based  on  11  Januaries,  1972-1982) 


f 

Daytime 
(12  L) 

Nighttime 
(00  L) 

0 

0.  311 

0.  380 

0.  1 

0.  016 

0.  015 

0.  2 

0.  026 

0.  013 

0.  3 

0.  026 

0.  022 

0.  4 

0.  030 

0.  018 

0.  5 

0.  025 

0.  021 

0.  8 

0.  033 

0.  024 

0.  7 

0.  038 

0.  027 

0.  8 

0.  058 

0.  037 

0.  9 

0.  137 

0.  099 

1.  0 

0.  300 

0.  345 

Mean 

0.  555 

0.  525 

-.umnia  rr/.i's  .aesc  estimates  of  mean  cloud  cover  for  all  45  layers,  for  the 
tion  iCi-ntfal  Europe,  northern  plains),  season  (July)  and  time  of  day 

■  I , 

-..iia-  kind  of  information  as  that  found  in  Table  5  was  entered  into  a  chart 
1  1  ind  cui'M.'s  drawn  las  shown).  The  curve  for  layers  beginning  at 
‘  -  fxt  r.i  poll  ted  downward  to  that  level,  thus  providing  an  estimate  of  the 

■I  ai  ,  i)\.-r  :it  that  level  (0.45).  Likewise,  other  estimates  were  made  for 
:  !■  •.•■1.-,,  and  the  envelope  drawn  for  the  levels.  Figure  13  thus  pro\  ides 
I  ;d  "'.er  it  all  levels  and  layers,  and  is  one  of  14  such  charts  that  have 
I  '  ti,  the  E  r Ac  data.  (The  probability  estimates  of  mean  cloud  cover  in 


•  iSed  in  the  examples  of  Figures  8-11.) 


WEAN  CuOUO  CO.EBS 


‘LOuD  COVER 

,M  ii.  ^^ERHa^JALYClS 

Fl  R  centb.m  EUROPE  . 
J&NuAR . ,  OOOCi 


Figure  13.  Chart  of  Mean  Cloud  Cover 
in  Central  Europe,  Southern  .Mountains, 
.lanuary  OOOOZ.  The  left-hand 
envelope  gives  the  mean  cloud  cover 
at  a  level  as  read  on  the  ordinate. 

Each  branching  curve  gives  the  mean 
cloud  clover  in  the  layer  between  the 
height  as  read  on  the  ordinate  and 
the  height  as  marked  on  the  curve 


southern  MTs 


Previously,  Gringoi'ten  had  proposed  that  the  mean  cloud  cover,  P^(h), 
within  a  layer  of  thickness  (h),  be  related  to  the  subtotals  of  mean  cloud  cover; 


P  (ID  and  P  (11+h),  bv  the  relation: 
o  o 


P  (hi  =: 
o 


P  (IPh)  -  P  (11) 
o  o 


If  this  is  a(  cepted,  the  problem  is  reduced  to  finding  a  suitable  expression,  or 
e.ytjressions,  foi’  tlie  parameter  C.  In  the  previous  study,  using  the  deBary  and 
jMoller  data  only,  C  could  be  made  a  function  of  h  and  H: 

CKIi.lli  .  -l.fifiF  i.!i4  yiT  +  0. 475h,  G,  h,  H  in  km 

F  545  111.1  yir  0.  475h,  G,  h,  H  in  feet  .  (26) 

'[’here  was  no  recognizable  relation  of  G  with  season  or  time  of  day. 

U  ith  the  later  E  fAC  data  the  relationship  between  the  parameter  G,  the  height 
of  the  (  loads  (II.  11*  h)  and  the  thickness  (h)  was  more  readily  given  as  follows: 


G(h)  F  HPO  +  0.  :i0h  ft 


-  6600  +  0.  66h  ft 


6000  +  1 . Olh  ft 


for  H  ^  2500 


2500  <  11  *  0500 


II  p  0500  . 


11.  Gringorten,  1.  I.  (1082)  (  liiiuitic  Probabilities  of  the  \  E  rtical  Distribution 
of  (  loud  Cover,  At- (jI.-'rK-82-007H,  .Al)  A11H753. 
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:is  m  tlu'  previousi  study,  the  value  of  G  could  not  be  systematically  /-elated 
til  tiiia-  Ilf  \cai-  ui-  time  of  day,  nor  could  location  be  clearly  associated  with  the 
.  line  of  (i.  riif  above  relations  are  considered  as  ustible  universally.  In  km: 


(!i) 

0.  27  1 

+  0.  3 Oh  km 

for  11  7- 

1  km 

2.  01 

a  0.  66 h  km 

1  < 

11  V  3  km 

1.  83 

a  1 . 0  Ih  km 

11  '• 

3  km  . 

(28) 

A  all  this  forniula  for  G.  liq.  (25)  should  provide  an  estimate  of  the  mean  cloud 

covi  i-  m  an\  lasiT,  for  ivhich  no  better  estimate  is  available. 

[o  see  if  Ikis.  i27I  and  (28)  are  suitable,  let  us  return  to  the  data  that  ETAC 

halpriuuled,  ttirouf^h  use  ifihe'I'DI-  1 4  (su  r  face  )  data,  as  exem  pli  fied  in  Table  3 .  The 

me  1/1  eloud  cover,  as  plotti-d  in  l  i^uire  12,  is  f/iven  only  for  tlie  whole  layer  from 

the  m  ound  to  levels  aloft.  The  application  of  Eq.  (27)  gives  the  mean  cloud  cover 

in  1  avers,  I'  (h).  Table  (j  is  a  record  of  these  estimates,  but  they  are  further 
o 

modified  to  reduce  the  estimate  in  a  layer  so  that  it  is  no  greater  than  the  estimate 
in  1  thicker  laver  that  contains  the  thinner  layer. 

U  hen  the  estimates  of  Table  6  are  entered  into  a  chart  (Figure  14)  it  appears 
that  there  a  l  e  Ihri’e  modal  heights  in  the  distribution  of  mean  sky  cover,  at  1000  ft, 
12,000  ft  and  30,000  ft  (low  stratus,  middle  clouds  and  cirrus?). 

Since  the  estimates  of  Table  6  or  1-  igure  14  are  not  verifiable  with  the  data  on 
hand,  tliev  can  be  accepted  only  if  they  are  satisfyingly  reasonable.  The  figures  in 
Table  5  Cil)  .Neph  estimates)  are  available  for  comparison  with  the  estimates  in 
T  dile  7.  Tile  KMSl)  of  the  estimates  of  mean  cloud  cover,  in  36  layers  between 
150  ft  an  1  55,000  ft.  was  found  to  he  0.066. 


Figure  14.  Model  Estimates  of 
Mean  C  loud  Cover,  for  Levels 
and  l.ayers,  for  Fort  Worth, 
Texas,  January,  0300L 


CLOUD  COVER 


6.  CONCLUSIONS 


The  immediate  goal  of  simulating  a  field  of  weather  conditions,  speiTfically 
cloud  cover,  has  been  achieved  through  use  of  the  "Three-Dimensional  Hoehni 
Sawtooth  Wave"  Model  (HD-BSWl.  It  is  possible  to  stochastically  simulate  a  field 
of  clouds,  in  the  horizontal  plane  or  in  a  vertical  cross-section.  The  simulation  is 
compatible  with  the  basic  statistics  of  cloud  cover,  including  the  climatic  mean  sky- 
cover  for  a  given  season  and  time  of  day,  frequency  distribution  of  cloud  covei- 
from  clear  to  overcast,  and  spatial  correlation.  This  r-eport,  however,  assumed 
that  the  parameters  of  the  3D-BSVV  model,  the  mean  cloud  cover  (B  )  and  the  scale 
distance  or  wavelength  {.\)  are  known.  Thi.s  report  mu.st  bo  followed  bv  another 
development  to  present  algorithms  for  the  determination  of  the  parameter  (A), 
using  published  climatological  information. 

The  basic  model  of  this  report  has  been  supplemented  with  supporting  models, 
to  extend  its  use  to  layers,  so  that  we  are  able  to  pi-ovide  simulation,  not  onlv  for 
the  total  sky  cover,  but  also  tor  a  specified  layer  in  the  atmosphere,  'i  (“t  tliei-e 
are  many  aspects  of  the  problem  that  must  await  development. 

Certainly  the  cloud  shapes  are  not  simulated  in  all  their  picturesque  variations, 
A  promising  idea  along  these  lines  is  the  simulation  of  the  "fractal"  sti-ucdure  of 
clouds.  The  characteristics  of  the  separation,  or  cloud-free  intervals,  should  be 
simulated,  and  tested.  Instead  of  a  2-dimensional  view,  a  single  line-of-ti-avid,  in 
and  out  <;f  clouds,  could  be  simulated  stochastically,  providing  a  realistic  time  and 
duration  of  exposure  of  a  flying  vehicle. 

■Still  to  be  accomplished  is  the  simulation  of  a  time-lapse  sequence  of  cloud- 
cover  events  in  three  dimensions.  Several  alternative  approaches  are  possible, 
including  the  modeling  of  the  Ornstein-Uhlenbeck  process.  An  investigation  into 
this  phenomenon  might  include  the  effects  of  terrain  and  other  geopliysical  leatures, 
\V  hen  time  changes  are  introduced,  the  simulation  of  advection,  and  cltjud  types, 
might  also  be  tried. 
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